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Abstract 

This paper studies a model of dark matter called wave dark matter (also known as 
scalar field dark matter and boson stars) which has recently also been motivated by a new 
geometric perspective by Bray [9]. Wave dark matter describes dark matter as a scalar field 
which satisfies the Einstein-Klein-Gordon equations. These equations rely on a fundamental 
constant T (also known as the "mass term" of the Klein-Gordon equation). In this work, we 
compare the wave dark matter model to observations to obtain a working value of T. 

Specifically, we compare the mass profiles of spherically symmetric static states of wave 
dark matter to the Burkert mass profiles that have been shown by Salucci et al. [3l] to 
predict well the velocity dispersion profiles of the eight classical dwarf spheroidal galaxies. 
We show that a reasonable working value for the fundamental constant in the wave dark 
matter model is T = 50 yr _1 . We also show that under precise assumptions the value of T 
can be bounded above by 1000 yr _1 . 



1 Introduction 

Ever since the first postulation of dark matter in the 1930's by Zwicky |40| , much evidence 
for the existence of dark matter has accumulated including the unexpected behavior in the 
rotation curves of spiral galaxies pp] , the velocity dispersion profiles of dwarf spheroidal galaxies 



31,37-39 , and gravitational lensing [13] . These and other observations support the idea that 
most of the matter in the universe is not baryonic, but is, in fact, some form of exotic dark matter 
and that almost all astronomical objects from the galactic scale and up contain a significant 
amount of this dark matter. Describing this dark matter is currently one of the biggest open 
problems in astrophysics [5|-[7] [T7][27]|30]|36] . 

In the last two decades, there has been substantial progress on describing the distribution 



of dark matter on the galactic scale. Navarro, Frenk, and White's popular model 26 , which 



resulted from detailed iV-body simulations, has been shown to agree well with observations 



outside the centers of galaxies 18, 38 . However, the Navarro-Frenk- White dark matter energy 
density profile also exhibits an infinite cusp at the origin, while observations favor a bounded 
value of the dark matter energy density at the centers of galaxies [14]|15]. This has prompted 



many astrophysicists to employ a cored profile, such as the Burkert profile [IT], which "cores 
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out the infinite cusp, to model the dark matter energy density in a galaxy. Resolving this 
"core-cusp problem" remains an important open problem in the study of dark matter. 

One possible solution to this problem, and potentially other astrophysical problems, is the 
introduction of a new model for dark matter. Recently, Bray has geometrically motivated the 
study of a scalar field satisfying the Einstein-Klein-Gordon equations as a viable dark matter 
candidate via constructing axioms for general relativity [9|j. The idea of using a scalar field 
to describe dark matter is not new. In fact, such a model has been seriously considered as a 
candidate for dark matter for more than two decades and has been shown to be in agreement 
with many cosmological observations 0[§|§[§[^[T9}[22}|25j[32}|35]. In most of these settings, 
these scalar fields are considered from a quantum mechanical motivation for the same Einstein- 
Klein-Gordon equations and go by the name scalar field dark matter or boson stars. However, 
due to the fact that the Klein-Gordon equation is a wave-type partial differential equation, we 
prefer the name wave dark matter. In this paper, we begin to test wave dark matter against 
observations at the galactic level. In particular, we seek a working estimate of the fundamental 
constant in the wave dark matter model, T, to be used in future comparisons to data. To do 
so, we will compare the simplest model defined by wave dark matter to models of dark matter 
that are already known to fit observations well. 

Salucci et al. recently used the Burkert profile to model the dark matter energy density 
profiles of the eight classical dwarf spheroidal galaxies orbiting the Milky Way. They found 
excellent agreement between the observed velocity dispersion profiles of these galaxies and those 
velocity dispersion profiles predicted by the Burkert profile 31 . This can be seen in Figure [TJ 
which we have reproduced exactly as it appears in the paper by Salucci et al. In what follows, 
we will show that a value of T = 50 yr _1 produces wave dark matter mass models that are 
qualitatively similar to the Burkert mass models found by Salucci et al. We will also show that 
under precise assumptions, comparisons to these Burkert profiles can be used to bound the value 
of T above by 1000 yr~ x . 



2 Burkert Mass Profiles 

The Burkert energy density profile models the energy density of a spherically symmetric dark 
matter halo using the function 

/W = 7 — : w 2 i 2\ 1 

(r + r c )(r z + r£) 

where po is the central density and r c is the core radius. Integrating this function over the ball 
of radius r with respect to the standard spherical volume form yields the Burkert mass profile 
as follows. 



M B {r)= / p B (s)dV R 3 

JB r (0) 
rr 

= 4-7T / s 2 [ib(s) ds 
Jo 
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Figure 1 : Observed velocity dispersion profiles of the eight classical dwarf spheroidal galaxies are 
denoted by the points on each plot with its associated error bars. The solid lines overlayed on 
these profiles are the best fit velocity dispersion profiles predicted by the Burkert mass profile. 
This figure is directly reproduced from the paper by Salucci et al. [31] and the reader is referred 
to their paper for a complete description of how these models were computed. 
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Burkert Mass Profile 




Figure 2: Plot of a Burkert mass profile. The inflection point is marked with an x. 
4W , gW J „ ds 



(s + r c )0 2 + r 2 ) 

Mb if) = 2 Wc 3 (in {Z±^j + \ In (^f) ~ arctan (2) 

A generic plot of a Burkert mass profile, Mg(r), defined to be the dark matter mass in the 
ball of radius r, is shown in Figure [2] We make a few remarks about the behavior of this mass 
function. 

Note that the behavior of the graph changes concavity at the inflection point r = ri p , which 
we have marked on the plot in Figure [2] with an x . Recalling from equation ^ the fact that 
Ms(r) is the integral over the interval [0, r] of the function 47rr 2 //s(r), we can compute this 
inflection point as follows. 

„3^2 



M B (r) = 4vrr /i B (r) = -— (3 

(r + r c ){r A + r£) 



Differentiating again yields 



-47r / 9 rg(r 4 -r 2 rg-2rrg) 
B{) (r + r c ) 2 (r 2 +r 2 ) 2 ' [> 



which has two complex zeros and two real zeros. The two real zeros are r = and 

/ 3+ (27 + 3 v / 78) 2/3 \ 



3(27 + 3^78) V3 



1.52r cs (5) 



the latter being the inflection point of the mass model. 

For r> Tip, the plot grows logarithmically due to the fact that the arctan term in equation (|2j) 
approaches a constant value as r — > oo. To describe the behavior when r <C rj p , we note that 
the Taylor expansion of Ms(r) centered at r = is as follows, 

M B (r) = ^p r 3 + O(r i ). (6) 
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Figure 3: Left: Plot of the Burkert mass profile for the Carina galaxy found by Salucci et al. 
|31| along with a mass plot of a wave dark matter static ground state, the cubic function which 
is the leading term of the Taylor expansion of the Burkert mass profile, and the quadratic power 

function — B ^ c ^ r 2 where r c is the core radius of the Carina galaxy. The x marks the location 

r c 

of the inflection point of the Burkert mass profile, while the vertical line denotes the location 
of the outermost data point for the Carina galaxy and is presented for reference purposes only. 
Right: Closeup of the plot on the left over the v interval [0,?"2p]. 



Thus for r <C Vi p , M#(r) is dominated by an r 3 term making the initial behavior cubic. 

In fact, several other models for dark matter mass profiles have similar initial behavior to 
the Burkert profile including a quadratic mass profile (which is not physical and is only included 
for the sake of comparisons) and wave dark matter mass profiles. In Figure [3j we have collected 
several mass models that have similar behavior inside r = ri v to the Burkert mass profile 



computed by Salucci et al. for the Carina galaxy 31 . While these models have similar behavior 
inside r = rj p , they are very different outside r = rj p . 

We have computed the inflection points of each of the Burkert mass profiles computed by 
Salucci et al. for the eight classical dwarf spheroidal galaxies [3JJ and have marked these points 
on a plot of each Burkert mass profile in Figure [4} We have constrained the viewing window of 
each plot to the range of data points collected. That is, we plot the Burkert mass profiles on 
the interval [0,r; as t], where ri ast is the radius of the outermost data point given by Walker et 



al. 38 39 1 for the observed velocity dispersion profiles. We have presented them in order from 
greatest to least according to the ratio of ri ast /ri p . 

In Table [TJ we have collected the defining parameters, po and r c , computed by Salucci et al. 
for the Burkert mass profiles which best predict the velocity dispersion profiles of each galaxy 
We have also collected the outermost data point, ri as t, of these velocity dispersion profiles 



39 , as well as our computations of the inflection point, rj p , and the ratio ri as t/ri p for each of 



the classical dwarf spheroidal galaxies. All quantities have been converted to geometrized units 
(the universal gravitational constant and the speed of light set to one) of (light)years for mass, 
length, and time. 
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Figure 4: Plots of the Burkert mass profiles computed by Salucci et al. of the eight classical 
dwarf spheroidal galaxies within the range of observable data. The inflection point is marked 
on each plot by an x . Carina and Draco have no inflection point marked because the inflection 
point for their Burkert mass profiles occurs outside the range of observable data. 
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Galaxy Name 
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ftp 


(yr) 


Tlast/ Tip 


Sextans 
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1.53 
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10 2 


3.26 
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2.32 
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14.05 


Leo II 
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1.88 
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1.37 
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8.57 
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10 3 
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1.84 
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3.01 
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9.19 
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3.03 
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10 3 


1.40 
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10 3 


2.17 
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1.10 
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1.16 
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10 3 


3.59 
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1.76 
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10 3 
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Ursa Minor 


1.83 
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8.01 
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10 2 


2.41 
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10 3 


1.22 
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10 3 


1.98 


Carina 


2.90 


X 


io- 


-16 


1.97 


X 


10 3 


2.84 


X 


10 3 


2.99 


X 


10 3 


0.95 


Draco 


8.19 


X 


io- 


-16 


2.11 
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10 3 


3.00 
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10 3 


3.20 


X 


10 3 


0.94 



Table 1 : Burkert mass profile data for the eight classical dwarf spheroidal galaxies converted to 
units of years for mass, length, and time. The parameters po and r c are those found by Salucci 
et al. for the best fit Burkert profiles [M], and ri as t is the radius of the outermost data point 
given by Walker et al. 38, 39 . Also included is the value of the inflection point, r?p, of the 
Burkert mass profile for each galaxy and the ratio of ri as t to ri p . 

3 Static States of Wave Dark Matter 

Now that we have presented mass profiles which model actual data well, we need to describe 
the wave dark matter models we will use to make our comparison. In the following, we present 
only the basic background information required to understand the model we use and refer the 
reader to [9j[l0] for more discussion on its motivation and successes thus far. 

Let (N,g) be a spacetime whose metric has signature ( — h ++). Let / : N — > C be a 
smooth complex-valued scalar field defined on the spacetime. Finally, let / and g satisfy the 
Einstein-Klein-Gordon equations, 

n a ( df®df + df®df (\df\\, f , 2 \\ (7 , 
G = 8vr/x l ^2 I ^2" + 1/1 J 9] (7a) 

□ s / = T 2 / (7b) 

where O g is the Laplacian with respect to the metric g. The parameter T is a fundamental 
constant of this system and its value must be determined in order to use these equations to 
model dark matter in the universe. On the other hand, the parameter /xo is not fundamental to 
the system and can be completely absorbed into / if desired. 

In [9], wave dark matter may be modeled as a real scalar field or as a complex scalar field. 
For the Einstein-Klein-Gordon system, a complex scalar field is equivalent to two real scalar 
fields. Complex scalar fields are more convenient because they have static spacetime solutions. 
Analogous "nearly static" real solutions can be achieved by redefining / for a static spacetime 
solution as v2 Re(/) whose spacetime metrics are very similar except for a high frequency 
oscillating pressure term which averages out to zero. Preliminary estimates suggest that this 
oscillating pressure effect might not be significant enough to make it measurable on physically 
relevant time scales. Hence, there may not be a testable difference between the predictions of 
real and complex scalar field dark matter at this time, though more thought on this question is 
well deserved. In either case, we will call these scalar field models "wave dark matter." 
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Since we wish to compare solutions of ([7]) to the spherically symmetric Burkert mass profiles 
computed by Salucci et al. 31 , we choose to work in spherical symmetry. In a recent paper 29 , 
Parry surveyed the well-known form of the metric of a general spherically symmetric spacetime 
in polar-areal coordinates, namely, 



(8) 



for real valued functions V and M and where da 2 = d6 2 + sin 2 6 dip 2 is the standard metric 
on the unit sphere. This metric has the following useful properties. The function M(t, r) is 
the Hawking mass of the metric sphere of radius r and time t. Under the Einstein equation, 
G = 8ttT, M is also the flat volume integral of the energy density term in the stress energy 
tensor. This motivates interpreting the function M(t, r) as the mass inside the metric sphere of 
radius r at time t. Finally, given the Einstein equation, in the low field limit, V is approximately 
the gravitational potential of the system. We refer the reader to [29] for detailed proofs of these 
facts. 

It is also shown in |29| that in spherical symmetry and using the metric ([8]), solving the 
Einstein-Klein-Gordon system Q reduces to solving the system 

(9a) 

(9b) 

(9c) 
(9d) 

To solve this system, we need boundary conditions. At the central value, we require all of 
the functions to be smooth. Since all of the functions are spherically symmetric, this implies 
that M r , V r , f r , and p r all vanish at r = for all t. We will also require that the spacetime 
be asymptotically Schwarzschild, that is, it approaches a Schwarzschild metric as r —> oo. 
Specifically, this implies that 

2V 4 2 (- 2M 




e 2V — >• « I 1 — — as r -> oo (10) 

D gs f ->■ T 2 / and / ->■ as r -> oo (11) 

where k > and gs is the appropriate Schwarzschild metric. Since / — > as r — > oo, M 
approaches a constant value m, which is the total mass of the system. 

Note that these boundary conditions ensure that as r — > oo, the metric g in equation Q 
becomes the Schwarzschild metric 

g s = -k 2 (l - ~) dt 2 + (l - ~) _1 dr 2 + r 2 da 2 . (12) 
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Thus k represents a scaling of the t coordinate in the standard Schwarzschild metric. The effect 
of this on our discussion is that V — > hiK as r — > oo. We have to numerically solve these 
equations and so in practice, we will impose these boundary conditions at an artificial right 
hand boundary point r max and solve the system on the r-interval [0,r ma:r ]. 

One of the simplest solutions to this system are those where the scalar field is of the form 



F(r) 



(13) 



f(t,r) 

where w 6 1 is a constant and F is real valued. 

Note that for / of this form, solving equation (11) for large r and requiring the solution to 
decay to yields that, for large r, F must satisfy 



+ 



1 



F ps 0. 



(14) 



F'+ WT 2 

\ V fC r J 

Requiring this condition on our system ensures that / appropriately decays to as r — > oo |28 



Solutions of the form in equation ( 13 ) produce static metrics and, once substituted into the 



system ^ , yield the following set of ODEs [28] , 



M' = 47rrVo 

V 
F' 
H' 



i + ^L e -2v 



! 2M \ \H\ 



I y 2 




Y2 



I -PI 



2M\ \H\ 



r J T 2 



H 



1 



2M\ 
r J 



-i 



lav 



with boundary conditions 
F(0) = 1, 



H(0) = 0, 



M(0) = 0, 



V(0) = V , 



(15a) 

(15b) 
(15c) 
(15d) 

(16) 



F'(r max ) + I Wr 2 -^ + 



V(r max ) - ^ In ( 1 - 2 Mi r 
2 V r, 




max 
max 



(17) 
(18) 



by equations (10) and (11). For simplicity, we set k = 1, which corresponds to the assumption 



on our choice of t coordinate that V goes to zero at infinity. A solution to these equations 
depends on the choice of the parameters T, fio, ui, and Vq. We solve a shooting problem for to 



and Vq to satisfy (17) and (18) leaving T and freely selectable. 



For each choice of T and /io, there are an infinite number of discrete finite mass solutions 
characterized by the number of zeros that F exhibits |28| . These are called static states. A 
static state with no zeros is called a ground state. With n zeros for n > 0, it is called an n th 
excited state. In Figure [5j we have plotted examples of F for a ground through third excited 
state. In Figure [6j we have presented the plots of the mass, M, corresponding to the plots of F 
in Figure [5j 
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Plot of F(r) for a Ground State Plot of F(r) for a 1 st Excited State 




Figure 5: Plots of spherically symmetric static state scalar fields (specifically the function F(r) 
in (|13|)) in the ground state and first, second, and third excited states. Note the number of 
nodes (zeros) of each function. 
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Figure 6: Mass profiles for a static ground state and first, second, and third excited states of 
wave dark matter. 
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3.1 Working Value of T 

As stated before, our goal is to find a value of T that is compatible with the Burkert mass 



profiles that Salucci et al. 31 computed to model the dark matter in the eight classical dwarf 
spheroidal galaxies. For T = 50 yr -1 , there is at least one wave dark matter static state that 
matches the Burkert mass profiles reasonably well. We have plotted such matches in Figure [7| 
Thus we have chosen to use 

T = 50 y^ 1 (19) 

as a working value of T in our future work with wave dark matter until we have the capability 
to make a more accurate approximation or precise measurement of this value. 

While we have chosen T = 50 yr -1 as a working value of T since it corresponds to wave 
dark matter models compatible with other well fitting models, the above does not constitute a 
precise measurement of the value of T. In the remainder of this paper, we will show that, under 
precise assumptions, the value of T can be bounded above, a first step towards obtaining a 
direct measurement of T. To do this, we first describe a few extra properties of these spherically 
symmetric static states. 

3.2 Hyperbolas of Constant T 



Parry showed in 28 how the parameters which define an n th excited state, T, and Vb, 

as well as the values of the total mass, m, and the radius, r^, called the half mass radius, for 
which M(rh) = m/2, are related to each other. The approximations which follow only apply to 
solutions in the long wavelength, low field limit, which is what is relevant for modeling galaxies. 
In particular, for constant T, there is a one parameter family of solutions for each n th excited 
state defined by the equations 

w»(T jA io) « Texp (cf requency ^j (20) 



^(^'^o) ~ Cpotential ~ yT" (21) 

m™(T,/, )«C^ ss T- 3 /Vo /4 (22) 

rMr,»o)*c? adius r-v 2 [i- 1/4 . (23) 

where the constants C" depend upon which state we wish to consider (i.e. they depend on 



n). Thus for constant T, the different possible excited states satisfying (16), (17), and (18) 
are defined entirely by the value of (j,q, which corresponds to the value of the energy density 
function, fi, as defined in |29j at the origin via the equation 



2 



Mo) = M 1 + Y2 e ~ 2V " ) • ( 24 ) 



The equations relevant to our discussion here are ([22]) and ( 23 ) . We have collected the values 
of C™ lflss and C^ adius from these two equations for the ground through fifth excited states as well 
as for the tenth and twentieth excited states in Table [2] A complete table of the values of all 
the constants CI 1 for these excited states can be found in |28|. 
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Figure 7: Static state mass profiles for T = 50 which are each a best fit to the Burkert profiles 
of the corresponding dwarf spheroidal galaxy. For T = 50, we picked an n th excited state whose 
best fit profile matched the Burkert profile qualitatively well. This shows that T = 50 is a 
reasonable working value of T. However, it does not imply that the actual value of T is 50 or 
that these galaxies are correctly modeled by the presented n th excited state. The units on T are 
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n 


^mass 


radius 





4.567 ±0.05 


0.8462 ± 0.004 


1 


10.22 ±0.10 


2.2894 ±0.009 


2 


15.81 ±0.16 


3.8253 ±0.014 


3 


21.37 ±0.22 


5.3994 ±0.018 


4 


26.91 ±0.27 


6.9860 ±0.022 


5 


32.42 ±0.33 


8.5606 ± 0.026 


10 


60.32 ± 1.18 


15.1357 ±0.039 


20 


116.62 ±2.57 


29.6822 ± 0.107 



Table 2: Values of the constants in equations (22) and (23) for the ground through fifth excited 



states as well as the tenth and twentieth excited states. We have given these values error ranges 
which encompass the interval we observed in our experiments, but it is possible that values 
outside our ranges here could be observed. However, we do not expect them to be outside by 
much if the discretization of r used in solving the ODEs is sufficiently fine. Note also that our 
values have less precision as we increase n. This is because as n increases, it becomes more 
difficult to compute the states with as much precision. 



As explained in (28), equations (22) and (23) imply that the product of m and rh does not 



depend on the value of /io, but only on T. Specifically, 

CmassC radius /nr\ 

mr h = ^2 ' ( 25 ) 

where we have suppressed the notation of n. If T is constant, then, because both C mass and 
Cradius are positive, the right hand side of this equation is some positive constant, k, and we 
have 

mrh = k (26) 

which defines a hyperbola. Thus, for a given n th excited state, all of the possible mass profiles 
for a constant value of T lie along a hyperbola. We illustrate this phenomenon in Figure [8j 

3.3 Fitting Burkert Mass Profiles 

With these properties of static state mass profiles in mind, we turn our attention to finding 
static state mass profiles that best fit the Burkert mass profiles computed by Salucci et al. [31| 
Given a Burkert mass profile, Mb, a value for T, and a specific state (i.e. value for n), we define 
for our purposes the best fit wave dark matter static state mass profile, M\y, as the one which 
minimizes the L 2 norm of the difference between these profiles, E, given by 

E = \\M B - M w \\\ 2 = [ '^(Mb - M w fdr. (27) 

Jo 

Of course, since we compute the static states numerically, we have to approximate this norm by 
an appropriate Riemann sum defined on a discretization of the interval [0,7*; as ^]. 
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Ground State Mass Profile 




Figure 8: Left: Plot of the mass profile of a ground state with its corresponding hyperbola of 
constant T overlayed. Any ground state mass profile that keeps the presented relationship with 
this hyperbola corresponds to the same value of T. Right: Examples of different ground state 
mass profiles corresponding to the same value of T. The corresponding hyperbola of constant T 
is overlayed. Notice that all three mass profiles have the same relationship with the hyperbola. 



To find this minimum, we first note that since T is fixed, we can write the total mass m and 



the value of /j,q in terms of a choice of via equations (25) and (23) respectively. Thus, we 
parameterize the different mass profiles of constant T, and hence E by r^, that is, E = E{rh). 
Furthermore, since all of the static state mass profiles of constant T lie on a hyperbola, there 
will be a value of Th that yields the minimum of E{rh). 

To make computing the best fits more uniform from galaxy to galaxy, we make the choice 
r/j = br c , where b > and r c is the core radius of the Burkert profile we wish to match, and vary 
the free parameter b. To compute which value of b produces a minimum of value of E(rh), we 
create a grid of values around an initial choice of b of the form [(b — step)r c , br c , (b + step)r c ] 
for some step > 0. Next we compute E(rh) for each of the values of Th and shift the grid, if 
necessary, so that it is centered on the value which yielded the smallest value of E(rh). If the 
grid shifts, we recompute E(rh) on the new grid and continue to shift, if necessary. Once the 
minimum E^r^) value occurs at the center of the grid, we keep that point as the center, but cut 
the step size in half. We then run this shifting procedure again for this smaller grid until the 
minimum is at the center and then we shrink again. We continue to shrink the step size until 
we get to a predetermined terminal value. We generally would run the procedure until the step 
size was less than or equal to 2~ 10 . 

We have already seen generic examples of the results of this matching procedure as the plots 
in Figure [7] are best fits for their respective galaxies, given an n th excited state, and the value 
T = 50 yr _1 . However, this best fitting procedure also provides a method of finding values of 
T, for T sufficiently large, which produce untenable matches to the Burkert profiles of Salucci 
et al. 1 31] This is the topic of the next section. 



3.4 Upper Bound for T 

To find an upper bound for T, we first need to explain how the static states change as T gets 



large. Equation (|25j) implies that for a given n th excited state, as T increases, the product mr^ 
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Ground State Mass Profiles for Leo II , l0 » Ground State Mass Profiles for Leo II 




Figure 9: Left: Ground state mass profiles of various values of T that are best fits to the Burkert 
mass profile found by Salucci et al. [3l] for the Leo II galaxy. The corresponding hyperbolas of 
constant T on which these profiles lie are also plotted. Ground states and their corresponding 
hyperbolas are drawn in the same color. Right: The same plots as in the left frame, but with 
the constant function which best fits the Burkert profile also plotted. Note that the best fit 
mass profiles approach this constant mass profile as T increases. 



decreases. The hyperbolas corresponding to smaller values of mr/, are those that lie closer to 
the mass and radius axes. 

Now consider the n th excited state mass profile that is the best fit to a Burkert mass profile 
for a given T. As T increases, the hyperbola to which this static state mass profile corresponds 
will get closer to the mass and radius axes, but since the mass profile must also minimize E(rh), 
the value of its mass will not tend to 0. Since mr^ tends to zero as T — ¥ oo, it must be instead 
that rh — >• as T — > oo. This implies that, as T increases, more of the constant portion of 
the best fitting n th excited state mass profile will be compared to the Burkert profile. Thus, as 
T — > oo, the best fitting n th excited state mass profile will limit to the constant function of r 
that best fits the Burkert profile under the same fitting criteria used for the static states. We 
illustrate this phenomenon in Figure [9j 

As T — > oo, the initial increasing region (i.e. the region before the constant portion) of the n th 
excited state mass profile that best fits a Burkert mass profile becomes more compressed. This 
initial region is where all of the dark matter mass is located. Thus as T increases, the dark matter 
corresponding to the best fit n th excited state extends out to smaller radii. However, observations 



suggests that dwarf spheroidal galaxies are dark matter dominated at all observable radii 20 
Thus the best fit n th excited state mass profiles for large T do not represent observations well 
and can be rejected. The question then is exactly when should we reject them. 

Since every static state has the initial increasing region just described, the best fit n th excited 
state for any value of T will be a better fit than the best fit constant function. Moreover, since 
this initial region becomes more compressed as T — > oo, for large T, the value of E(rh) for the 
best fit n th excited state mass profile increases monotonically as T — > oo approaching the value 
of E for the best fitting constant function. 

This suggests a criteria for when to reject values of T. We will reject a best fit n th excited 
state mass profile, and hence its corresponding value of T, as an untenable model of the dark 
matter mass if its value of E{rh) is greater than some prescribed fraction of the value of E 
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Figure 10: The Burkert mass profile found by Salucci et al. [3l] for the Sextans galaxy. The best 
fit static state mass profiles for a ground through fifth excited state, tenth excited state, and 
twentieth excited state all lying on the same hyperbola are overlayed on the plot. The hyperbola 
here satisfies the rejection criteria for all of the different static states represented in the plot, 
thus all of these static states correspond to an upper bound on the value of T for Sextans for 
their respective value of n (i.e. the set of n th excited states). Note how close together all of the 
states are. This is due to the fact that the majority of their profiles which are being compared 
to the Burkert mass profile is the common and constant portion of the profiles. 



for the best fitting constant function. We chose to use 80%. Explicitly, we use the following 
rejection criteria. 

Rejection Criteria 3.1. Given T, n, and a Burkert mass profile Mb, let M\y be the spherically 
symmetric n th excited state mass profile corresponding to T that best fits Mb, that is, that 



minimizes E from equation (21) along the hyperbola defined by the value of T and equation 
(25). Let Eyy be the value of E for this mass profile. Furthermore, let Mq be the constant 
function which best fits Mb, also by minimizing the corresponding function E, and let Eq be 
the value of E for the constant function. Reject the given value T as a tenable value for this 
fundamental constant if 

Ew — SEq- 

In other words, any fit that is less than 20% better than the best fitting constant function 
of r is rejected as a bad fit. 

For each of the eight classical dwarf spheroidal galaxies and n £ {0,1,2,3,4,5,10,20}, we 
computed values of T that yielded n th excited state mass profiles that best fit that galaxy's 
Burkert profile which were rejected by the above criteria. All of the values of T above those 
computed are also rejected because they produce mass profiles even closer to the constant 
function. In Table [3j we have collected these upper bounds of T. In Figure 10, for the galaxy 



Sextans, we present best fit static state mass profiles for the ground through fifth, tenth, and 
twentieth excited states for which E(rh) is more than 80% of the value of E for the best fit 
constant function. 
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Galaxy \ State 











1 






2 






3 




Sextans 


T 


< 


160 


T 


< 


394 


T 


< 


633 


T 


< 


875 


Leo II 


T 


< 


234 


T 


< 


576 


T 


< 


926 


T 


< 


1279 


Fornax 


T 


< 


35 


T 


< 


87 


T 


< 


139 


T 


< 


192 


Leo I 


T 


< 


57 


T 


< 


141 


T 


< 


226 


T 


< 


312 


Sculptor 


T 


< 


49 


T 


< 


121 


T 


< 


194 


T 


< 


268 


Ursa Minor 


T 


< 


82 


T 


< 


202 


T 


< 


325 


T 


< 


449 


Carina 


T 


< 


84 


T 


< 


207 


T 


< 


333 


T 


< 


459 


Draco 


T 


< 


45 


T 


< 


111 


T 


< 


179 


T 


< 


246 



Galaxy \ State 




4 




5 




10 




20 


Sextans 


T 


< 


1116 


T 


< 


1356 


T 


< 


2460 


T 


< 


4789 


Leo II 


T 


< 


1632 


T 


< 


1983 


T 


< 


3597 


T 


< 


7003 


Fornax 


T 


< 


245 


T 


< 


297 


T 


< 


538 


T 


< 


1048 


Leo I 


T 


< 


398 


T 


< 


484 


T 


< 


877 


T 


< 


1706 


Sculptor 


T 


< 


342 


T 


< 


416 


T 


< 


754 


T 


< 


1467 


Ursa Minor 


T 


< 


572 


T 


< 


695 


T 


< 


1261 


T 


< 


2455 


Carina 


T 


< 


586 


T 


< 


712 


T 


< 


1292 


T 


< 


2514 


Draco 


T 


< 


314 


T 


< 


382 


T 


< 


692 


T 


< 


1347 



Table 3: Upper bound values for T corresponding to poor best fits of the Burkert mass profiles 
for each of the classic dwarf spheroidal galaxies. The values in each column for each galaxy 
should be interpreted as an upper bound on the value of T, under the approximations explained 
in the paper, if that galaxy is best modeled by an n th excited state. The units on T are yr _1 . 
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We observe from Table [3] that the upper bound values of T increase as we increase the state 
we consider. This is due to the following. The rejected values of T correspond to rejected 
hyperbolas of constant T, and hence constant mr^. Furthermore, the only qualitative difference 
between any two n th excited state mass profiles is the number of ripples in the initial increasing 
region of the profile. For large T, the majority of a best fit n th excited state mass profile that 
is compared to the Burkert profile is the constant region which is shared by static state mass 
profiles for any n. Thus the hyperbola corresponding to a rejected best fit ground state is close 
to the hyperbola corresponding to a rejected best fit n th excited state for any n. In particular, 
there is a hyperbola of constant rar/,, for which the corresponding best fit n th excited state mass 
profiles for any n are rejected by the above criteria. Then, since C^ ass and C^ adius appear to 
monotonically increase as n increases (see Table [2j , by equation ( 25 ) , we would expect the same 



behavior for the value of T in order for mr^ to remain constant, which is what we observe in 
Table H 

Thus if a dwarf spheroidal galaxy is correctly modeled by a twentieth excited state or less, 
then an overall upper bound on the value of T would be the upper bound corresponding to the 
twentieth excited state. The least upper bound corresponding to the twentieth excited state 
over all eight galaxies is that value for the Fornax galaxy, which yields approximately that 

T < 1000 yr -1 . (28) 

4 Utilized Approximations 

Now that we have presented our results, we list here the important approximations made in this 
paper which led to these results and explain briefly why we make them. 

Approximation 1: Dark matter is correctly described by the wave dark matter model. 

Approximation 2: Dark matter halos around dwarf spheroidal galaxies are spherically sym- 
metric. 

Approximation 3: The dwarf spheroidal galaxies used in this paper are in a state of dynamical 
equilibrium. 

Approximation 4: The spacetime metrics describing these dwarf spheroidal galaxies are static. 

Approximation 5: Wave dark matter predicts outcomes qualitatively similar to those of spher- 
ically symmetric static state solutions to the Einstein-Klein-Gordon equations. 

Approximation 6: The Burkert mass profiles computed by Salucci et al. |31| fit the observa- 
tional data very well. 

Approximation 7: The spacetime is in the low field limit, that is, M <C r. 

Approximation 8: The spacetime is asymptotically Schwarzschild. 

Approximation 1 is used because we are testing the wave dark matter model against obser- 
vations. Approximation 2 is a common approximation for dwarf spheroidal galaxies and is also 
necessary because we are comparing the wave dark matter model to the spherically symmet- 
ric Burkert mass profile. Approximation 3 seems to be consistent with observations of dwarf 
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spheroidal galaxies at least out to large radii |12 31 . Approximation 6 is reasonable given 



Figure [TJ Approximations 7 and 8 are standard when modeling galaxies. 

Approximations 4 and 5 are used to simplify the types of solutions to the Einstein-Klein- 
Gordon equations we consider. We note here, however, that there is a question of the stability 
of the spherically symmetric static state solutions. It is known that, if the corresponding total 
mass is not too large, the ground state is stable under perturbations (21.33 but that, on their 
own, the excited states are not [I] regardless of their mass. However, it has also been shown that 
a coupling of an excited state with a ground state can produce a stable configuration |3|. We 
hypothesize that luminous matter distributions coupled with combinations of static states will 
produce a stabilizing effect allowing for more dynamically interesting systems to be physically 
plausible. 



5 Conclusions 

To summarize the results of this paper, we have drawn effectively two conclusions, which we list 
here. 

Conclusion 5.1. Given Approximations 1 through 8, a value of T which yields one or more 
spherically symmetric static state mass profiles which match well the best fit Burkert mass profiles 
computed by Salucci et al. [31] for each of the eight classical dwarf spheroidal galaxies is 

T = 50 yf 1 . 

Conclusion 5.2. Given Approximations 1 through 8 and Rejection Criteria \3.1\ if the dark 
matter halos of all of the eight classical dwarf spheroidal galaxies are correctly modeled by 20 th 
excited states or less, then 

T < 1000 yr' 1 . 

For the interested reader, the Matlab code used for this paper to generate the spherically sym- 
metric static states and to compute the best fits to a Burkert profile can be found on Bray's Wave 
Dark Matter Web Page at http : //www. math. duke . edu/~bray/darkmatter/darkmatter .html, 
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